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Abstract 

The  diffusion-kinetics  model  for  the  inters  tween  interleukin-2  and 

each  of  its  T-cell  surface  receptors  (IL-2o:  and  IL-2#)  is  presented.  This  model  is 
unique  in  that  it  considers  both  three  dimensional  ligand-receptor  interactions 
and  two  dimensional  interactions  between  cell  surface-bound  species.  Elemen¬ 
tary  rate  laws  are  developed  for  initial  encounters,  rebounding  interactions,  and 
dissociations  of  free  ligands  and  receptors  according  to  the  method  of  Waite  [1] 
and  Waite  and  Stewart[2].  Analogous  rate  laws  are  written  for  membrane  bound 
species  which  undergo  similar  initial  associations,  rebounding  interactions,  and 
dissociations.  A  set  of  kinetic  equations  is  proposed  for  a  system  consisting  of 
two  independent  monovalent  receptors  and  one  monovalent  ligand,  simulating 
the  interaction  of  the  IL-2ce  and  IL-2/J  receptors  of  the  human  T-cell  with  the 
lymphokine  interleukin-2.  Autocrine  and  paracrine  growth  and  combinations 
of  the  two  are  studied  by  modifying  the  appropriate  experimental  parameters. 
Experimental  associative  and  dissociative  rate  constants  are  determined  for 
important  T-cell  surface  species, 
keywords:  Interleukin-2,  diffusion,  signalling,  kinetics 
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1  Background 


The  cellular  membrane  is  the  barrier  of  the  cell  at  which  many  biologically  important 
interactions  take  place.  These  interactions  include  transport  of  materials  in  and  out 
of  the  cell  and  the  transduction  of  signals  to  the  cell  [3].  This  signal  transduction 
instructs  the  cell  to  replicate,  synthesize  a  compound,  etc. 

Under  the  currently  accepted  fluid  mosaic  model  [4],  the  cellular  membrane  con¬ 
sists  of  a  phospholipid  bilayer  with  the  hydrophilic  heads  pointing  outward  and  the 
hydrophobic  tails  on  the  inside,  making  a  waterproof  barrier.  Many  proteins  can  be 
associated  with  this  lipid  bilayer  through  various  weak  bonding  interactions.  Recep¬ 
tors  are  proteins  that  are  in  effect  solvated  in  the  cellular  membrane  and  held  in  place 
by  weak  hydrophilic  interactions,  whereby  nonpolar  parts  of  the  protein  molecule  in¬ 
teract  with  the  hydrophobic  portion  of  the  cellular  membrane.  The  actual  receptor 
site  of  the  protein  molecule  is  a  region  with  a  special  shape  that  exists  on  the  outside 
of  the  cellular  membrane.  The  entire  receptor  molecule,  being  solvated  in  the  lipid 
bilayer,  is  free  to  diffuse  in  two  dimensions  and  randomly  does  so,  in  search  of  the 
specific  substrate  for  that  site.  The  substrate  may  be  a  ligand  diffusing  towards  the 
cell  surface  in  the  extra-cellular  medium,  or  it  may  be  bound  to  another  receptor 
diffusing  within  the  two-dimensional  membrane. 

The  interaction  of  T-cells  and  the  lymphokine  interleukin-2  (IL-2)  is  one  such 


receptor-substrate  interaction.  The  T-cells  are  key  participants  in  the  cell-mediated 
immune  response  of  the  body,  and  their  action  is  mediated  by  interleukin-2. 

When  a  T-cell  is  presented  an  antigen  and  recognizes  it  through  a  signal  trans¬ 
duced  through  the  cellular  membrane,  it  begins  to  produce  the  lymphokine  interleukin- 
2,  which  is  a  molecular  communicator  and  growth  factor  for  various  cells  of  the  im¬ 
mune  system  [5j.  Interleukin-2  is  necessary  for  the  continued  growth  of  the  T-helper 
cell  and  for  the  growth  of  the  T-cytotoxic  cell  which  lyses  cells  containing  the  antigen 
(6,  7]. 

The  growth  of  T-cells  can  be  initiated  by  interleukin-2  in  one  of  two  fashions.  In 
autocrine  growth,  interleukin-2  produced  by  a  particular  T-cell  becomes  bound  to 
the  IL-2  receptors  on  that  same  cell,  causing  it  to  continue  its  growth  and  production 
of  interleukin-2.  The  interleukin-2  produced  by  a  particular  T-cell  may  also  escape 
this  cell’s  surface  and  enters  the  bulk  solution.  From  there  it  may  bind  to  the  IL-2 
receptors  of  another  T-cell,  causing  it  to  grow  and  produce  interleukin-2.  This  is 
known  as  paracrine  growth  [8]. 

In  actual  biological  systems  T-cell  growth  most  likely  occurs  by  some  combination 
of  these  two  methods.  Both  the  autocrine  and  paracrine  growth  mechanisms  are  of 
interest  here  because  they  increase  the  number  of  cells  which  can  recognize  and  defend 
against  a  specific  antigen. 

The  nature  of  intracellular  signaling  is  not  of  interest  to  this  study.  Only  the 
kinetics  of  association  and  dissociation  of  interleukin-2  and  its  T-cell  receptors  that 


lead  to  the  signal  transduction  are  significant  to  this  study.  The  molecular  signal  that 
is  generated  once  the  interleukin-2-receptor  complex  is  formed  is  also  not  of  interest 
here. 

It  is  currently  believed  that  the  IL-2  receptor  is  not  merely  a  single  protein  but 
rather  two  proteins  that  must  both  be  bound  to  interleukin-2  in  order  to  stimulate  the 
continued  production  of  interleukin-2  and  the  growth  of  T-cells  [9].  The  interleuldn- 
2  receptor  is  believed  to  consist  of  two  proteins  (10],  the  receptor  (IL-2a),  which 
is  a  protein  with  a  molecular  weight  of  55kD,  and  the  receptor  (IL-2/?),  which  is  a 
protein  with  a  molecular  weight  of  75kD.  It  is  interesting  to  note  that  these  receptors 
themselves  do  not  appear  until  the  T-cell  has  been  activated  by  an  antigen  becoming 
bound  to  the  T-cell  [11].  It  is  known  [12,  13]  that  IL-2q  is  a  low-affinity  receptor  for 
interleukin-2,  having  a  K*  of  lxlO“8M  and  that  IL-2/?  is  a  medium-affinity  receptor 
for  interleukin-2,  having  a  Kd  of  lxlO~9M.  It  is  also  known  that  interleukin-2  rapidly 
associates  and  dissociates  with  IL-2a  but  associates  and  dissociates  slowly  with  IL-2/? 
[10]. 

The  actual  signal  transduction  occurs  when  interleukin-2  is  bound  to  both  the 
IL-2a  receptor  and  the  IL-20  receptor,  forming  a  ligand-receptor  complex  [5].  It  is 
interesting  to  note,  however,  that  interleukin-2  binds  rapidly  to  the  complex  and  yet 
dissociates  slowly  from  the  complex  [10]. 

The  objective  of  the  current  study  is  to  quantitatively  study  the  proposed  mech¬ 
anisms  for  receptor-ligand  signaling  in  the  IL-2  problem.  There  is  an  ongoing  debate 
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over  which  mechanism  is  the  correct  model  for  the  formation  of  the  interleukin-2-T- 
cell  surface  receptors  complex. 

The  affinity  conversion  mechanism  for  the  interaction  of  T-  cells  and  interleukin-2 
is  explained  as  follows.  Interleukin-2  produced  by  an  activated  T-cell  that  is  wan¬ 
dering  in  solution  encounters  an  IL-2  receptor  that  is  moving  randomly  in  the  two- 
dimensional  cellular  membrane.  If  the  IL-2  has  correct  orientation  and  energy,  all 
of  which  can  be  accounted  for  in  a  single  association  probability  parameter,  it  will 
rapidly  bind  to  the  IL-2  receptor.  This  portion  of  the  complex  continues  to  be  an¬ 
chored  in  the  cellular  membrane  by  IL-2  and  will  move  randomly  within  the  mem¬ 
brane.  If  the  interleukin-2/receptor  complex,  while  diffusing  randomly  in  the  cell 
membrane,  encounters  the  opposite  type  of  IL-2  receptor  with  sufficient  energy  and 
proper  orientation  then  the  complete  interleukin-2  -  receptor  complex  is  formed  and 
an  intracellular  signal  is  transmitted,  instructing  the  T-cell  to  begin  producing  more 
IL-2  and  increase  its  rate  of  growth.  The  affinity  of  the  IL-2a  receptor  is  effectively 
converted  to  that  of  the  IL-2$  receptor  when  the  complex  is  formed  in  this  manner. 

The  preformed  heterodimer  mechanism  produces  the  signaling  complex  in  a  dif¬ 
ferent  manner  [14,  15,  16,  17,  18].  Before  interleukin-2  becomes  involved  in  the 
process,  the  IL-2a  receptors  and  IL-2^  receptors  must  first  find  each  other  in  the 
two-dimensional  cell  membrane  and  form  a  complex.  This  complex  has  been  termed 
the  high-affinity  interleukin-2  receptor,  having  a  K d  of  lxlO-nM  [17].  Once  the  IL- 
2a/IL-2/3  heterodimer  has  been  formed,  interleukin-2  must  diffuse  through  the  bulk 


solution  and  bind  to  the  heterodimer  in  a  ligand-surface  association. 

The  area  of  controversy  in  regards  to  this  mechanism  is  the  question  of  whether  the 
receptors  are  truly  separate  in  the  cellular  membrane  or  if  a  heterodimer  of  both  the 
IL-2  and  the  IL-2  exists  as  the  receptor  [14, 15,  16,  17,  18].  This  recent  experimental 
work  purports  to  have  shown  that  a  heterodimeric  receptor  for  IL-2  is  the  proper 
model,  although  the  evidence  presented  there  is  still  not  conclusive.  The  problem  with 
these  studies  is  that  they  attempt  to  make  use  of  equilibrium  techniques  to  describe 
what  is  indeed  a  kinetics  problem.  It  is  the  intent  of  this  study  to  use  theoretical 
computer-modeling  of  kinetics  to  determine  which  of  these  two  proposed  mechanisms 
is  the  one  which  best  describes  the  formation  of  the  interleukin-2  -receptor  complex. 

2  Introduction 

2.1  General  Description  of  the  Model  and  Study 

The  study  of  the  kinetics  of  biological  systems  has  become  a  useful  tool  for  under¬ 
standing  the  actual  interactions  involved  in  such  systems  [14].  The  present  study 
uses  theoretical  kinetics  to  model  the  formation  of  IL2-IL2a  complexes  and  IL2-IL2/? 
complexes  on  the  surface  of  the  human  T-cell. 

Recently  much  work  has  been  done  on  the  interleukin-2  system  [14,  15,  16,  17, 
18].  Of  particular  interest  is  the  question  of  the  presence  of  preformed  receptor 
heterodimers.  Much  of  this  work  has  supported  the  existence  of  such  preformed 


heterodimers  but  detailed  theoretical  kinetic  results  have  not  yet  been  provided. 

It  is  the  purpose  of  the  present  study  to  use  the  method  of  Waite  [1]  to  model  the 
kinetics  of  the  interleukin-2  system.  The  diffusion-kinetics  interactions  of  ligands  with 
cell-membrane  bound  receptors,  including  the  formation  of  ligand- receptor  complexes 
which  lead  to  signal  transduction,  are  described  by  a  set  of  differential  equations  which 
are  obtained  by  combining  the  appropriate  rate  terms  for  each  of  the  species  in  the 
model.  There  are  analogous  rate  terms  for  free-ligand  cell  surface  interactions  and 
membrane-bound  species  interactions  with  parameters  for  each  adjusted  accordingly. 
The  proposed  species  for  this  study  are  listed  in  Table  1. 

This  present  study  involves  both  a  simple  and  a  complex  model  of  the  diffusion- 
kinetics  of  the  interleukin-2  system.  The  simple  model  is  a  preliminary  step  to  the 
eventual  theoretical  diffusion-kinetics  solution  to  the  interleukin-2  preformed  het¬ 
erodimer  problem.  This  model  considers  the  interaction  between  interleukin- 2  and 
the  IL-2a  receptor  and  the  interaction  between  interleukin-2  and  the  IL-2/9  receptor 
but  not  the  interactions  between  IL-2a  receptors  and  IL-2/?  receptors. 

The  complex  model  is  used  to  actually  solve  the  interleukin-2  preformed  het¬ 
erodimer  problem.  It  allows  the  IL-2a  receptors  and  IL-20  receptors  to  interact  and 
possibly  form  the  heterodimer  IL-  2a/IL~20  which  can  then  either  dissociate  or  bind 
interleukin-2  to  form  the  signaling  complex. 

Ligands  bound  non-specifically  to  the  cell  surface  (glycocalyx)  are  included  in 
both  models.  Parameters  such  as  diffusion  coefficients,  probabilities  of  binding,  and 


lifetimes  of  involved  diffusing  species  are  arrived  at  theoretically  [2]  or  by  fitting 
theoretical  data  empirically  to  actual  experimental  data. 
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2.2  Quantitative  Aspects  of  the  System 

The  system  used  in  this  study  consists  of  a  volume  of  solution,  V,  containing  n,*// 
cells  of  radius  Rceii-  Each  cell  surface  has  n-sites  of  radius  r  with  a  diffusion  coefficient 
D  [1].  For  the  simple  model  these  sites  consist  of  species  2-11  as  listed  in  Table  1A. 
For  the  complex  model  these  sites  consist  of  species  2-17  as  listed  in  Tables  1A  and 
IB. 

Each  cell  can  be  assigned  varying  initial  numbers  of  species  and  varying  probability 
factor,  diffusion  coefficients,  and  dissociation  constants.  The  concentrations  of  all 
species  involved  in  the  system  are  more  conveniently  written  in  terms  of  number 
densities  rather  than  molar  concentrations  as  all  calculations  are  based  on  the  number 
of  species  per  cell. 

The  surface  number  density  (m~2)  of  a  particular  species,  /?,  is  simply  the  number, 
rip,  of  that  particular  species  present  divided  by  the  surface  area  of  the  cell  [1] 

N0  =  n0/AvR7eell  (1) 

In  this  model  all  of  the  &  species  are  sites  where  specific  binding  can  take  place. 
These  sites  include  species  2-17  as  defined  previously  in  Table  1. 

The  fraction  of  the  cell  surface  occupied  by  these  sites  is  [1] 


fff  -  NpTTp  (2) 

As  shown  by  Waite  [1],  the  non-specific  binding  sites  on  the  glycocalyx  of  the  cell 
do  not  have  number  densities.  The  fraction  of  cell  surface  available  for  non-specific 
binding  is  given  by 


=  (3) 

B 

where  the  summation  is  over  all  of  the  possible  binding  sites,  /?,  defined  above  [1]. 

Bulk  free  ligand  particles  are  also  present  in  this  system.  They  are  described  in  a 
similar  manner  as  the  surface  bound  species  in  that  they  have  radius  r;  and  diffusion 
coefficient  t>i  [1].  The  volume  number  density  of  ligands  (m“3)  is 

N,  =  M,Na  •  1000  (4) 

where  M/  is  the  molar  concentration  of  ligand  and  NA  is  Avogadro’s  number  [1]. 

2.3  Dimensional  Aspects  of  the  Diffusion-Kinetics  Model 

Diffusional  rate  constants  are  based  on  classical  kinetics  with  some  important  modifi¬ 
cations.  These  differences  are  necessitated  by  the  fact  that  the  receptors  and  ligands 
involved  in  the  theoretical  kinetics  model  do  not  undergo  strictly  three  dimensional 


interactions. 
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An  encounter  of  two  bodies  in  solution  can  be  modeled  by  the  following  second 
order  rate  law  [1]: 


rate  =  Jb[A][f?]  (5) 

where  k  is  the  rate  constant  for  the  reaction  in  m3s-1  and  [A]  and  [£]  are  the  con¬ 
centrations  of  species  A  and  B  in  m-3,  respectively. 

A  problem  arises  when  this  type  of  rate  law  is  applied  to  a  ligand-receptor  inter¬ 
action.  This  is  because  the  ligand  is  in  solution  and  has  three  degrees  of  freedom,  but 
the  receptor  is  solvated  in  the  cellular  membrane  and  has  only  two  degrees  of  freedom. 
For  the  normal  solution  case,  the  second-order  rate  constant  would  be  expected  to 
have  units  of  m3s-1  as  shown  above. 

Waite  [1]  has  shown  that  if  the  rate  of  disappearance  of  a  substrate  is  being 
measured,  then  the  rate  constant  must  be  expressed  in  terms  of  m7s~l  [10].  Similarly 
if  the  rate  of  disappearance  of  a  membrane  bound  species  is  being  measured,  then 
the  rate  constant  must  be  expressed  in  terms  of  m3s_1  [1].  It  is  convenient  then 
to  define  a  conversion  between  these  rate  constants  for  surface  species  and  those  for 
solution-phase  species.  Waite  [1]  has  defined  this  conversion  as 


k3o  =  Ck7D  (6) 

where  C  =  to  R^n^n /V ,  which  is  the  ratio  of  the  total  surface  area  of  all  the  cells 
in  the  system  to  the  total  volume  of  the  system. 
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3  Derivation  of  Ligand- Surface  Kinetic  Terms 

3.1  Ligand- Surface  Associations 


The  total  flux  J  of  initial  encounters  of  free  ligand  with  a  cell  surface  of  radius  R^u 
is  given  by  [1] 


Jl-eell  =  ixDlRcellNt 


(7) 


The  total  rate  of  disappearance  of  free  ligand  due  to  encounter  with  the  cell  surface 
is  [1] 

rate  =  W  (8) 

where 

Kl‘£i  =  A™eeilD,Rcell/V  (9) 

Upon  initial  encounter  of  free  ligand,  l,  and  a  cell  receptor  site,  /?,  one  of  three 
events  is  possible.  The  ligand  can  bind  to  the  receptor  site  with  probability  Pp-i,  the 
ligand  can  escape  to  the  bulk  medium  and  become  free  once  again  with  probability 
(1  —  Pfi-i)  Pete,  or  the  ligand  can  rebound  from  the  receptor  site  and  enter  a  non-free 
state  with  probability  (1  —  Pp-i)  (1  —  Petc)-  These  interactions  are  represented  by  the 
following  kinetic  steps  [1]: 
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1-0 


/  +  M 


1---0  > 

1  +  0 


(10) 


The  term  PM c  is  the  probability  that  a  ligand  can  diffuse  from  the  cell  surface 
into  the  bulk  without  subsequent  re-encounters  with  the  cell  surface  and  is  equal  to 
Ti/ {r edi  +  r/)  [1].  Each  of  these  probability  factors  modifies  the  rate  constant  for  a 
given  type  of  interaction. 

The  explicit  term  of  the  total  rate  of  appearance  of  species  resulting  from  bulk 
to  surface  associations  is  the  total  flux  of  free  ligand  at  the  cell  multiplied  by  the 
fraction  of  the  surface  area  occupied  by  receptor  site  [1]: 


rate  =  J  fp/V  =  i^riccuDiRceurpNiNp/V  (11) 

The  modification  of  this  rate  law  by  the  appropriate  probability  factors  gives  the  rate 
of  appearance  of  species  1—0  and  l  •  •  •  0.  No  rate  term  is  needed  for  the  case  where 
the  result  is  1  +  0  because  there  has  been  no  net  change  in  the  numbers  of  0  or  / 
present. 

It  would  appear  that  as  a  species  formed  it  could  itself  become  a  receptor  site 
and  eventually  lead  to  the  introduction  of  many  very  complex  species  into  the  model. 


The  simplifications  of  these  species  take  the  general  form 


(0  —  l)  —  l  sa  0  —  l  and  g  -  l 
(0  •  •  •  /)  —  l  0  —  1  and  g  •  •  •  / 

(12) 

(0  —  /) •  •  •  /  «  0  —  l  and  g  '•  -  I 
(0  **'/)••*/  ft*  0  •  •  ■  /  and  g  •  •  •  1 

3.2  Ligand-Surface  Dissociations 

All  bound  species  eventually  undergo  dissociation  by  first-order  kinetics  with  rate 
constants  designated  The  rate  law  for  such  dissociations  is 


rate  =  kp-i[Np-t)  (13) 

i 

Two  events  are  possible  upon  dissociation:  the  ligand  can  remain  associated  with 
the  cell  surface  through  rebounding  interactions  with  probability  1  —  Pete  or  the  ligand 
can  escape  to  the  bulk  medium  without  subsequent  re-encounters  with  the  cell  surface 
with  probability  Petc.  The  dissociation  of  bound  species  can  be  represented  by  the 
following  elementary  kinetic  steps  [1]: 


1-0 


(14) 


1---0 

1  +  0 

The  first-order  dissociation  rate  constant  of  a  bound  complex  is  related  to  the 


half-life  ,Ti/2,  of  the  bound  complex  as  shown  below 
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tor  depending  on  what  type  of  dissociative  event  is  occurring. 

3.3  Ligand-Surface  Rebounding  Interactions 

Localized  free  ligands  which  are  loosely  associated  with  the  cell  surface  are  all  possible 
species  0  ■  •  •  /.  These  species  undergo  re-encounters  with  the  cell  surface  with  a  total 
average  frequency  of  [1] 

vtot  —  2(1  —  In  2  )£>//*■?  (16) 

Re-encounters  at  the  surface  will  occur  either  at  the  site  of  origin,  0,  or  some  new 
location  on  the  cell  surface,  a.  For  simplicity,  rebounding  encounters  at  the  site  of 
origin  are  assumed  to  have  the  same  probability  parameters  as  if  the  ligand  had  been 
of  the  free  type  [1].  The  possible  outcomes  of  rebounding  associations  both  at  the  site 
of  origin,  0,  or  new  cell  surface  sites,  a ,  are  represented  by  the  following  elementary 
kinetic  steps  [1] 

0-1 
0  •  -  •  l 
0  +  1 


0---1 


(17) 
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0  +  a-l 

/9«--Z  +  or-»  «  0  +  a-->l  -  (18) 

0  +  a  + 1 

►  ^ 

The  simplifications  of  complex  species  described  above  are  invoked  here  to  avoid 
introducing  unnecessary  complex  species  into  the  model. 

The  total  re-encounter  frequency,  described  above,  includes  all  possible  re-encounter 
sites.  This  can  be  divided  into  the  frequency  of  re-encounter  at  the  site  of  origin,  vp, 
and  the  frequency  of  re-encounter  at  all  other  surface  sites,  i/re#  [1],  This  frequency 
of  re-encounter  at  the  site  of  origin  is  given  by 

VS  =  ^ S(n  +  l)!(n+l)!2=*+>  1  _eXp  (4^  +  l)r?)]  (19) 

and  the  frequency  of  re-encounters  beyond  the  site  of  origin  is  given  by 


Vrt,  =  Vtoi  ~  vp 


These  frequency  terms  become  the  rate  constants  for  these  rebounding  interactions. 


4  Derivation  of  Membrane- Bound  Diffusion-Kinetics 
Terms 

4.1  General 


In  general,  all  of  the  bulk  to  surface  terms  derived  above  can  be  applied  to  the  case 
of  membrane  bound  interactions.  In  this  case,  however,  initial  encounter  rates,  re- 


encounter  frequencies,  and  escape  probabilities  are  all  based  on  results  obtained  from 
the  two-dimensional  diffusion  equation  [1].  The  membrane  bound  steps  also  need 
only  consider  interactions  between  species  that  may  successfully  bind,  since  it  is  only 
these  interactions  which  lead  to  the  formation  of  new  species,  whereas  all  bulk  to 
surface  interactions  give  rise  to  new  species  [1]. 

4.2  Membrane- Bound  Associations 

Consider  the  free  surface  species  u>  and  7  , diffusing  randomly  in  the  two-dimensional 
matrix  of  the  cell  membrane.  Each  species  has  diffusion  coefficient,  Du  or  and  ra¬ 
dius,  rw  or  r^.  The  term  bulk  free  represents  species  that  are  not  part  of  an  interacting 
pair  and  are  beyond  the  boundary  of  interactions,  X  [1]. 

The  total  rate  of  membrane  bound  encounters  between  species  u  and  7  is  given 
by  [1] 


rate  =  k„^[Nu}[Ny] 


with  a  rate  constant  of  [1] 


=  2  ir(D^  +  Dw)/]n 


(21) 


(22) 


In  this  case,  L  is  simply  the  average  distance  between  particles  on  the  cell  surface  [1]. 
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These  encounters  can  be  represented  by  the  following  elementary  kinetic  steps  [1] 


u>  —  7 
u  •  •  •  7 
i*>  +  7 


(23) 


Analogous  to  the  ligand  surface  case,  the  formation  of  complex  u  —  7  is  governed 
by  the  probability  factor  the  formation  of  the  species  u  •  •  •  7  is  governed  by 

the  probability  factor  (1  —  PMe)  (1  —  Pu- r)i  and  the  return  of  both  particles  to  the 
bulk  after  encounter  is  represented  by  (1  —  Pu-i)  Pt»e  [1]»  In  these  terms  Pu-. y  is  the 
probability  of  binding  per  encounter  and  PMe  is  the  probability  of  escape  to  the  bulk 
state  and  is  given  by  [1] 

p“‘  =  ln2/‘"(^)  (24) 


As  in  the  ligand  surface  case,  a  newly  formed  bound  species  may  also  serve  as  a 
free  surface  species  u  or  7.  Simplifications  are  made  in  this  case,  analogous  to  those 
preformed  in  the  ligand  surface  case,  in  order  to  eliminate  highly  complex  species 
from  the  model. 


4.3  Membrane>Bound  Dissociations 

Complexes  formed  from  bound  surface  species  also  undergo  dissociation  by  first-order 
kinetics.  Each  surface  dissociation  occurs  according  to  the  rate  law  where  ku-y  is  the 
rate  constant  defined  by  Waite  [1].  The  rate  of  this  dissociation  is  given  by 
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rate  =  kd„toc,  [A^]  ( 25) 

where  kdutoc.  is  the  first-order  dissociation  rate  constant. 

Two  outcomes  are  possible  upon  surface  dissociation:  the  species  u>  and  7  can 
re-encounter  one  another  to  form  the  species  w  •  •  •  7  or  the  species  u  and  7  can 
recede  beyond  the  distance  L,  becoming  free  surface  species  without  undergoing  sub¬ 
sequent  associative  re-encounters.  Formation  of  the  species  u  •  •  •  7  is  governed  by 
the  probability  parameter  (1  —  Petc)  and  formation  of  free  surface  species  is  governed 
by  the  probability  parameter  Pete  [1].  The  elementary  kinetic  steps  representing  such 
dissociations  are  [1]: 


u>  —  7 


u)  •  •  • 

u>  +  7 


(26) 


4.4  Membrane-Bound  Rebounding  Associations 


The  total  frequency  of  re-encounter  for  surface  species  v  •  •  7  is  given  by  [1] 


^  =  °-460((^) 


(27) 


The  elementary  kinetic  steps  representing  these  rebounding  interactions  are  [1] 


u>  —  0 


u  —  0 

^  (aJ  •  •  •  P 

U)  +  P 


(28) 
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Rebounding  encounters  leading  to  bound  species  u>  —  7  are  governed  by  the  probabil¬ 
ity  term  Pu—n  rebounding  encounters  leading  to  the  formation  of  rebounding  species 
u> •  •  *7  are  governed  by  the  probability  parameter  (1  —  Pw-n)  (1  —  P^c),  and  rebound¬ 
ing  encounters  leading  to  the  production  of  free  surface  species  are  governed  by  the 
probability  parameter  (1  -  Pw_7)  PMe  [1]. 

5  Method  of  Writing  Rate  Laws 

Rate  laws  governing  the  time  rate  of  change  of  bulk  ligand  and  possible  surface  species 
were  written  using  the  appropriate  bulk  to  surface  kinetic  terms  and  membrane- 
diffusion  kinetic  terms.  These  terms  where  then  translated  into  FORTRAN  code. 
The  rate  law  for  any  given  species  is  formed  by  summing  all  of  the  rate  law  terms 
for  the  mechanistic  steps  that  contribute  to  the  appearance  or  disappearance  of  that 
particular  species. 

A  brief  example  of  this  is  provided  by  examining  rate  terms  governing  the  time 
rate  of  change  of  the  species  resulting  from  the  bulk  to  surface  association  of  a  ligand 
molecule  and  a  cell  surface  receptor  site,  /?. 

•  Mechanistic  Step  for  Ligand-Surface  Association  [1] 
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a)l-P 

l  +  0  -  ■  b)  l  •  •  •  0  (29) 

c)l  +  fi 

i  4 

It  is  of  interest  to  note  that  no  rate  terms  are  needed  for  step  c  because  there  is 
no  net  change  in  the  number  of  free  ligand  or  receptor  site 

•  Rate  Terms  for  Ligand-Surface  Association  [1] 

The  following  terms  govern  the  appearance  of  the  products  indicated  in  the  steps 
given  above. 


a)  rate  =  Pi-iJi-ctiifi/VC 

(30) 

b)  rate  =  (1  -  />-,)(  1  -  PeM-aufi/VC 
Identical,  negative  terms  are  written  for  the  disappearance  of  the  reactants,  except 

that  the  conversion  from  three  dimensions  to  two  dimensions  ( C  =  4R^eHneen/V)  is 

not  needed  in  this  case  for  the  disappearance  of  the  ligand. 

•  Rate  Terms  Translated  into  FORTRAN 

As  shown  above,  the  rate  term  for  this  case  is  the  total  flux  of  encounters  between 
bulk-ligand  and  a  cell  surface  (J)  multiplied  by  the  appropriate  probability  factor 
divided  by  the  total  volume  of  the  system  and  the  three  dimensional  to  two  dimen¬ 
sional  conversion  factor.  The  value  of  the  flux  was  used  frequently  in  the  FORTRAN 
code  and  hence  labelled  as  a  constant  variable,  rl. 


a)  FORTRAN  term  =  +rl*frac(2,y)*pil/c 

b)  FORTRAN  term  =  +rl*frac(2,y)*(l-pil)*(l-pesc)/c 

The  term  frac(2,y)  is  simply  the  fraction  of  the  cell  surface  covered  by  the  i 
receptor.  Corresponding  numbers  are  assigned  to  other  species  in  the  model  as  shown 
in  Table  1.  is  the  probability  of  association  for  the  bulk  to  surface  interaction  of 
free  ligand  and  an  i  receptor  and  Pete  is  the  probability  of  escape  from  the  cell  surface 
without  re-encounter.  Experimentally  important  parameters  are  defined  in  Table  2. 

The  resulting  FORTRAN  program  is  run  in  order  to  generate  the  kinetic  data. 
The  coupled  differential  equations  are  solved  using  the  DIVPAG  subroutine  which 
is  based  on  the  method  of  Gear.  The  number  of  differential  equations  needed  is 
equal  to  the  number  of  species  present  in  the  model,  eleven  for  the  simple  model  and 
seventeen  for  the  complex  model.  The  data  is  analyzed  using  the  ttro  spreadsheet 
package  QUATTRO,  version  4.0. 

6  Experimental  Kinetic  Runs  Leading  to  the  So¬ 
lution  of  the  Preformed  Heterodimer  Problem 

6.1  The  Simple  Model  with  Identical  Receptors 

Kinetic  runs  were  performed  using  the  eleven  species  given  in  Table  1A.  It  is  of 
interest  to  note  that  these  species  include  no  possible  interactions  between  the  two 
types  of  receptors.  This  is  to  insure  that  the  system  is  well  behaved  with  no  negative 


concentration  values  and  a  constant  sum  of  the  numbers  of  both  i  and  j  ligands 
between  all  species  where  these  ligands  are  involved. 


25 


Two  types  of  runs  were  carried  out  in  order  to  elucidate  the  behavior  of  two 
types  of  natural  intercellular  signaling:  autocrine  growth  and  paracrine  growth.  The 
autocrine  growth  runs  allowed  no  initial  free  bulk  ligand  and  100000  initial  non* 
specifically  bound  ligands  in  addition  to  50000  each  of  i  and  j  receptors.  The  ma¬ 
jority  of  the  runs  describe  paracrine  growth  because  this  is  the  major  mechanism  of 
intercellular  signaling.  These  runs  had  50000  each  of  i  and  j  receptors  and  a  bulk  free 
ligand  concentration  of  lxl0-9M. 

6.2  The  Complex  Model  with  Identical  Receptors 

i 

The  complexity  of  the  system  was  increased  with  the  addition  of  i  and  j  receptor 
interactions.  This  change  gave  rise  to  the  additional  species  presented  in  Table  IB. 
The  addition  of  these  species  was  important  in  showing  whether  or  not  the  system 
continued  to  be  well-behaved  as  the  types  of  interactions  were  made  increasingly 
complicated.  As  before,  ’well-  behaved’  is  defined  to  mean  no  negative  numbers  of 
species  and  a  constant  sum  of  i  and  j  receptors  regardless  of  their  form,  i.e.  complexed 
or  unbound,  bound  to  a  particular  receptor,  etc. 

Three  types  of  kinetic  runs  were  investigated  in  order  to  elucidate  the  behavior 
of  the  system  for  each  type  of  intercellular  signaling.  As  before  both  autocrine  and 
paracrine  growth  were  studied,  with  the  addition  of  a  realistic  autocrine  and  paracrine 


growth  system  in  which  free  bulk  ligand  and  non-specifically  bound  ligand  is  provided 
at  the  beginning  of  the  kinetic  run.  The  autocrine  runs  have  3000  non-specifically 
bound  ligand  and  1000  each  of  i  and  j  receptors.  The  paracrine  runs  have  an  initial 
ligand  concentration  of  1x10 ~9M  and  1000  each  of  i  and  j  receptors.  The  combination 
runs  have  initial  concentration  of  ligand  of  lxl0-9M,  1000  each  of  i  and  j  receptors, 
and  3000  non-specifically  bound  ligand. 

6.3  The  Complex  Model  with  Non- identical  Receptors 

The  actual  interleukin-2  system  is  modeled  using  differing  receptors.  The  IL-2a 
receptor,  represented  in  the  model  by  the  i  receptor,  has  a  higher  binding  affinity  for 
interleukin-2  than  the  IL-2/9  receptor,  represented  in  the  model  by  the  j  receptor. 

As  in  the  complex  model  with  identical  receptors,  autocrine  growth,  paracrine 
growth,  and  the  combination  of  the  two  are  modeled  using  the  same  concentrations 
of  bulk  ligand,  non-specifically  bound  ligand,  and  i  and  j  receptors  as  described  for  the 
case  of  the  complex  model  with  identical  receptors.  It  was  these  runs  which  provided 
the  theoretical  kinetics  solution  to  the  preformed  heterodimer  problem. 


7  Results  and  Discussion 


7.1  The  Simple  Model 

7.1.1  Default  Data 

The  first  step  in  understanding  the  complete  paracrine  growth  mechanism  was  to 
obtain  data  showing  the  time  rate  of  change  of  the  amounts  of  the  species  of  interest 
under  the  initial  conditions  in  the  default  data  set.  Representative  ligand-receptor 
association  curves  for  this  data  set  are  given  in  Figure  1.  This  run  provided  the 
baseline  data  for  comparison  to  other  sets  of  data  obtained  using  modified  parameters, 
which  will  ultimately  be  the  tool  used  to  elucidate  the  overall  mechanism. 

7.1.2  Change  of  Cell  Radius 

At  all  times  the  larger  cell  (Reeii  =  lxl0_5m)  has  a  greater  number  of  non-specifically 
bound  ligands,  due  to  the  fact  that  its  surface  area  is  100  times  larger  than  the  default 
cell  (Reeii  =  lxl0~6m)  (Figure  la).  Initially,  the  number  of  complexes  of  either  type 
is  greater  for  the  default  cell.  This  is  because  the  surface  density  of  receptors  is  higher 
in  the  default  cell,  making  it  easier  for  free  ligand  to  ’’find”  an  i  or  j  receptor.  Figure 
lc  shows  that  at  3  milliseconds,  the  larger  cell  begins  to  have  more  i-1  complexes 
than  the  default  cell.  This  is  counter  to  the  explanation  offered  above  due  to  the 
surface  area  change  and  therefore  another  mechanism  must  predominate:  membrane- 
bound  association.  There  is  more  non-specifically  bound  ligand  on  the  surface  of 
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the  larger  cell  due  to  its  larger  surface  area,  providing  a  greater  surface  density  of 
non-specifically  bound  ligand  in  the  larger  cell.  This  makes  it  easier  for  the  surface 
bound  i  and  j  receptors  to  associate  with  membrane  bound  ligand  in  the  larger  cell. 
The  point  in  time  where  the  default  cell  begins  to  lag  behind  the  larger  cell  in  the 
number  of  bound-receptor  ligand  complexes  is  the  point  in  time  where  the  mechanism 
switches  from  predominantly  bulk  to  surface  associations  to  predominantly  membrane 
diffusion  associations. 

The  predominance  of  the  membrane-bound  association  is  further  shown  by  com¬ 
parison  of  association  rate  constants  for  the  complex  i-1  for  cells  having  radii  of  1 
micron  and  10  microns  (Figure  lb).  The  rate  constants  are  determined  using  the 
method  of  initial  slopes  and  are  given  and  compared  to  values  obtained  by  Landgraf, 
et  al.,  in  Table  3  [14]. 

7.1.3  Elimination  of  Ligand-Surface  Associations  for  I-L 

The  only  way  to  form  the  bound  receptor-ligand  complex,  i-1,  was  through  membrane 
diffusion  associations.  The  ligand  was  allowed  to  form  bound  receptor-ligand  com¬ 
plexes  with  the  j  receptor  via  both  bulk  to  surface  association  and  cell  membrane 
association. 

The  ligand-receptor  curve  for  this  case  is  shown  in  Figure  2a.  The  difference 
between  the  numbers  of  i-1  in  the  case  where  ligand  surface  interactions  have  been 
eliminated  and  the  default  case  becomes  progressively  less  significant  with  time.  At 
all  times  the  number  of  i-1  formed  by  both  bulk  to  surface  association  and  membrane- 
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diffusion  association  is  greater  than  the  number  of  i-1  formed  by  membrane-diffusion 
association  alone.  At  longer  times,  the  number  of  i-1  complexes  formed  when  ligand- 
surface  associations  have  been  eliminated  approaches  the  number  of  i-1  complexes 
formed  in  the  default  case. 

It  is  of  interest  to  note  that  the  rate  of  formation  of  j-1  complexes  when  IL-2- 
receptor  ligand  surface  associations  have  been  removed  is  greater  than  the  rates  of 
formation  of  i-1  and  j-1  complexes  under  the  default  conditions.  The  rate  constants 
are  determined  by  the  method  of  initial  slopes  (Figures  2b, c)  The  effective  number 
of  j-1  complexes  formed  increases  in  this  case  because  the  j  receptor  does  not  have  to 
compete  with  the  i  receptor  for  free  ligand  in  bulk  to  surface  associations. 

i 

7.1.4  Elimination  of  Membrane-Bound  Associations  for  I-L 

In  this  case  the  only  way  to  form  the  bound  receptor-ligand  complex,  i-1,  is  via  a  bulk 
to  surface  association  mechanism.  As  in  the  trials  where  ligand-surface  interactions 
were  eliminated  for  the  receptor  and  interleukin-2,  the  ligand  was  allowed  to  complex 
with  the  j  receptor  via  both  mechanisms. 

The  differences  between  the  number  of  i-1  in  the  case  where  membrane-bound  as¬ 
sociations  have  been  eliminated  and  the  default  data  becomes  increasingly  significant 
with  time  (Figure  2a).  At  larger  times  this  difference  increases  to  the  point  where 
the  i-1  complexes  formed  in  the  default  case  vastly  outnumber  the  i-1  complexes  due 
to  the  fact  that  the  i-1  complex  can  be  formed  via  both  mechanisms  in  the  default 
case  and  the  i-1  complex  can  be  formed  by  the  bulk  to  surface  mechanism  alone  in 
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the  other  case. 

It  is  interesting  to  note  that  the  rate  of  formation  of  j-1  complexes  formed  when 
membrane-bound  associations  have  been  eliminated  is  greater  than  the  rate  of  for¬ 
mation  of  i-1  and  j-I  complexes  formed  under  default  conditions.  This  is  because  the 
j  receptor  does  not  have  to  compete  with  the  i  receptor  for  ligand  in  membrane- 
diffusion  associations,  which  is  not  allowed  experimentally  for  the  i  receptor.  The 
rate  constants  were  determined  by  the  method  of  initial  slopes  as  shown  in  Figures  2 
b  and  c. 

7.1.5  Elimination  of  Ligand-Surface  Associations  for  I-L  and  Change  of 
the  Cell  Radius 

As  in  the  case  when  Rati  =  lxl0“6m  and  ligand-surface  associations  have  been  elim¬ 
inated  (Figure  2a),  the  production  of  i- 1  complex  initially  lags  behind  the  production 
of  i-1  complex  in  the  default  case  (Figure  3a)  when  Rail  —  lxl0~5m  and  ligand-surface 
associations  have  been  eliminated.  At  longer  times  the  production  of  i-1  complex  for 
this  case  surpasses  that  of  the  default  case.  The  larger  cell  has  more  non-specifically 
bound  ligand  that  can  be  converted  to  bound- receptor  ligand  complex  by  membrane- 
bound  associations,  which  are  the  only  way  that  the  complex  i-1  may  form.  The 
association  rate  constants  for  this  case  are  calculated  using  the  method  of  initial 
slopes  Figure  3b  and  are  given  in  Table  3. 


7.1.6  Elimination  of  Membrane-Bound  Associations  and  Change  of  the 
Cell  Radius 

As  in  the  case  when  Rctii  —  lxlO-6m  and  membrane  bound  associations  have  been 
eliminated  (Figure  2a),  the  rate  of  production  of  i-1  complex  for  this  case  initially 
matches  the  rate  of  production  of  i-1  complex  for  the  default  case.  At  long  times, 
however,  the  production  of  i-1  complex  in  this  case  cannot  keep  pace  with  the  pro¬ 
duction  of  i-1  in  the  default  case  (Figure  3a). 

The  larger  cell  size  makes  it  harder  for  the  ligand  to  find  the  receptors  on  the  cell 
surface  since  the  overall  number  density  of  the  receptors  is  decreased  by  a  factor  of 
100.  Ligand-surface  associations  are  the  only  way  to  form  the  complex  i-1  and  thus 
the  fact  that  the  larger  cell  has  more  non-specifically  bound  ligand  is  immaterial 
Thus,  the  rate  of  i-1  complex  production  for  the  case  of  larger  cell  elimination  of 
membrane-bound  interactions  is  less  than  the  rate  of  i-1  complex  production  in  the 
default  case,  as  shown  by  the  rate  constants  in  Table  3.  These  rate  constants  were 
calculated  using  the  method  of  initial  slopes  as  shown  in  Figure  3c. 

7 .1.7  Relationship  Between  Non-specifically  Bound  Ligand  and  Receptor- 
Ligand  Complexes 

The  ligand  receptor  association  curves  for  these  trials  are  given  in  Figure  4.  When 
the  bulk  to  surface  mechanism  is  predominant,  the  numbers  of  non-specifically  bound 
ligands  and  ligand-receptor  complex  behave  independently.  At  longer  times  the  dif- 
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fusion  kinetics  mechanism  takes  predominance  and  the  number  of  non-specifically 
bound  ligands  stays  relatively  constant  as  the  number  of  bound  i  receptor-ligand 
continues  to  grow  because  the  non-specifically  bound  ligands  collide  with  unbound 
receptors  to  form  bound  receptor-ligand  complexes. 

As  the  number  of  bound  receptor- ligand  complexes  approaches  its  maximum  value 
(50003  for  these  trials)  the  numbe.  non-specifically  bound  ligands  increases  dramat¬ 
ically.  There  are  progressively  fewer  and  fewer  free  receptors  of  either  type  to  interact 
with  the  non-specifically  bound  ligand  because  most  of  them  have  been  used  up  in 
receptor-ligand  complexes.  The  only  place  that  a  ligand  can  bind  to  the  cell  is  the  cell 
membrane  itself.  The  result  is  the  rapid  increase  in  the  number  of  non-specifically 
bound  ligand-cell  membrane  complexes.  This  result  is  even  more  exaggerated  in  a 
system  containing  cells  with  radii  of  lxl0~sm.  Even  more  of  the  cell  surface  is  avail¬ 
able  for  the  non-  specific  binding  of  ligand  to  the  cell  membrane  after  all  of  the  free 
i  and  j  receptors  have  been  consumed. 

7.1.8  Elimination  of  Non-specifically  Bound  Ligand 

These  studies  were  carried  out  with  cells  of  radius  lxl0-6m  and  lxl0~5m  in  order  to 
elucidate  the  relationship  between  non-specifically  bound  ligand  and  receptor-ligand 
complexes  shown  in  Figure  4  and  discussed  above.  In  both  cases,  when  ligand-surface 
interactions  have  been  eliminated  (Figures  5  a,b)  the  number  of  i-1  and  j-1  receptor- 
ligand  complexes  reached  an  equilibrium  number  of  2.50  before  receding  back  to 
zero  because  of  dissociation.  The  initial  amount  of  non-specifically  bound  ligand 
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was  quickly  used  by  the  i  and  j  receptors  in  the  formation  of  i*l  and  j-1  complexes. 
Dissociation  of  these  complexes  gave  bulk  free  ligand  instead  of  non-specifically  bound 
ligand  because  the  probability  of  formation  of  non-specifically  bound  ligand  was  set 
to  0. 

When  membrane-bound  associations  were  eliminated  (Figures  5a, b)  the  number 
of  receptor-ligand  complexes  almost  reached  the  maximum  value  of  50003  for  these 
experimental  conditions.  The  complexes  i-1  and  j-1  are  produced  by  solely  by  the 
ligand-surface  association  mechanism.  Thus  non-specifically  bound  ligand  is  necessary 
for  the  formation  of  bound  receptor- ligand  complexes  and  has  no  effect  on  the  ligand- 
surface  mechanism. 

7.1.9  Half-Lives  of  I-L  and  G-L  Complexes 

The  half-life  of  the  i-1  bound  receptor  ligand  complex  was  obtained  by  setting  elimi¬ 
nating  both  ligand-surface  and  membrane-bound  associations  in  the  model.  Therefore 
once  a  complex  had  dissociated  there  was  no  way  for  re-association.  The  first  order 
rate  constants  of  dissociation  were  obtained  from  the  first  order  plot  shown  in  Figure 
6  for  non-specifically  bound  ligand  and  the  complex  i-1  for  cells  with  a  radius  of  1 
micron  and  10  microns.  These  rate  constants  and  half-lives  are  listed  in  Table  3  and 
compare  favorably  to  actual  experimental  work  [14]. 
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7.1.10  Modification  of  the  Diffusion  Coefficient  of  the  Ligand 

In  the  case  where  Di  =  lil0-5m2s-1,  the  diffusion  constant  of  the  ligand  is  in¬ 
creased  to  the  point  that  the  rate  of  formation  of  bound-receptor  ligand  complexes  is 
greater  than  that  seen  in  the  default  data  where  Di  =  lxl0'"1®ni2s"1.  The  membrane- 
diffusion  mechanism  predominates  and  produces  bound-receptor  ligand  complex  more 
rapidly  because  the  ligand  is  better  able  to  move  in  the  two-  dimensional  matrix  of 
the  cell  membrane  (Figure  7a). 

When  Di  =  lxlO-15m2s-1,  the  diffusion  constant  of  the  ligand  is  decreased  to 
the  point  that  it  moves  much  slower  while  solvated  in  the  cellular  membrane.  This 
makes  it  harder  for  receptors  and  non-  specifically  bound  ligand  to  associate.  The 
great  increase  in  rate  of  complex  production  seen  in  the  default  data  from  membrane- 
diffusion  associations  is  absent.  The  bulk  to  surface  mechanism  predominates  and  can¬ 
not  produce  large  numbers  of  bound  complexes  as  quickly  as  the  membrane-diffusion 
mechanism  (Figure  7a).  Changing  the  cell  radius  to  10  microns  gave  similar  results 
in  both  cases  as  shown  by  Figure  7b. 

7.1.11  Study  of  Autocrine  Growth  Via  Elimination  of  FVee  Ligand  and 
Ligand-Surface  Associations 

The  general  behavior  exhibited  when  autocrine  growth  is  the  growth  mechanism  for 
the  system  is  that  the  complex  i-1  forms  at  the  same  rate  whether  or  not  ligand- 
surface  associations  are  included  in  the  model.  This  is  shown  by  the  ligand-receptor 
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association  curves  given  in  Figure  8.  Tlie  insignificant  difference  between  the  number 
of  i-1  complex  in  the  default  autocrine  case  and  in  the  case  where  ligand-surface 
interactions  have  been  removed  from  the  model  shows  that  the  primary  mechanism 
for  autocrine  growth  is  membrane-diffusion  associations.  Bulk  to  surface  associations 
are  insignificant  in  autocrine  growth  because  little  interleukin-2  manages  to  escape 
the  cell  surface  and  become  the  bulk  free  ligand  necessary  to  give  rise  to  ligand-surface 
associations. 

7.2  The  Complex  Model  with  Identical  Receptors 

7.2.1  Experimental  Conditions 

All  runs  for  the  complex  model  with  identical  receptors  were  conducted  with  1000 
of  each  type  of  receptor  present  on  the  cell  surface.  The  autocrine  growth  runs  had 
3000  non-specifically  bound  ligand  present  on  the  cell  surface  and  there  was  no  bulk 
free  ligand  present  in  the  system.  The  paracrine  growth  runs  had  no  non-specifically 
bound  ligand  present  initially  and  a  bulk  free  ligand  concentration  of  1x10 ~9M.  The 
combined  autocrine-paracrine  growth  mechanism,  which  most  closely  simulates  the 
actual  interleukin-2  system,  had  3000  non-specifically  bound  ligand  present  on  the 
cell  surface  and  an  initial  bulk  free  ligand  concentration  of  lxl0"9M.  Each  receptor 
has  identical  probability  parameters  so  that  associations  for  each  take  place  at  the 


same  rate. 
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7.2.2  Autocrine  Growth  Case 

Figure  9  shows  the  ligand-receptor  association  curves  for  autocrine  growth.  In  the 
default  data  case  the  signaling  complex,  represented  by  (i*j)*l,  reaches  an  equilibrium 
value  of  320  complexes  at  a  time  of  10  milliseconds.  Other  important  species  include 
the  i-1  and  j-1  complexes,  which  reach  an  equilibrium  value  of  620  complexes  at  5 
milliseconds,  and  the  preformed  heterodimer  which  reaches  a  maximum  of  90  com¬ 
plexes  at  1  millisecond.  This  plot  shows  that  the  complex  i-1  is  formed  faster  than 
the  preformed  heterodimer  i-j. 

In  the  case  where  ligand-receptor  interactions  are  eliminated  ( P\—i  —  0,  Pj—i  — 
0,  Pn-i  =  0,  P7j-i  =  0)  a  completely  different  set  of  ligand-  receptor  curves  are  ob¬ 
tained.  The  complexes  i-1  and  j-1  are  not  present  during  the  course  of  the  run  because 
the  probabilities  of  association  either  by  ligand-surface  associations  or  membrane- 
bound  associations  are  0.  The  equilibrium  for  the  (i-j  )-l  signaling  complex  has  changed 
with  respect  to  the  default  data  case.  The  equilibrium  number  of  the  signaling  com¬ 
plex  achieved  is  much  greater  than  that  of  the  default  data  case  and  reaches  a  max¬ 
imum  of  990  complexes  and  an  average  equilibrium  value  of  975  complexes.  The 
equilibrium  value  is  achieved  at  0.5  seconds,  much  later  than  the  default  data  case. 

The  complex  i-j  in  this  case  reaches  a  greater  maximum  that  in  the  default  data 
case  because  there  are  no  i-1  or  j-1  complexes  being  formed.  Thus  all  of  the  i  and  j 
receptors  are  available  for  preformed  heterodimer  formation. 

The  case  when  preformed  heterodimer  formation  is  eliminated  (Pi-j  —  0)  shows 


behavior  that  is  similar  to  the  default  data  set.  The  signaling  complex  (i-j)-l  reaches 
an  equilibrium  value  of  320  complexes  at  a  time  of  10  milliseconds,  as  was  the  case 
for  the  default  data  set.  The  i-1  complex  ligand-receptor  association  curve  for  this 
case  also  matches  that  of  the  default  data  case. 

7.2.3  Autocrine  and  Paracrine  Growth  Case 

Figure  10  shows  the  ligand-receptor  association  curves  for  the  autocrine  and  paracrine 
growth  system.  This  most  closely  models  the  actual  behavior  of  the  immune  system 
where  cells  receive  signals  via  both  mechanisms.  Here  the  behavior  of  the  default 
case  closely  matches  that  of  the  autocrine  growth  only  system.  The  (i-  j)-l  signaling 
complex  reached  an  equilibrium  value  of  320  complexes  at  10  milliseconds  and  the  i-1 

i 

and  j-l  complexes  reached  an  equilibrium  value  of  620  complexes  at  9  milliseconds. 
As  was  the  case  with  the  autocrine  growth  only  runs,  only  a  small  number  of  the 
preformed  heterodimer  (i-j)  was  present  and  quickly  disappeared. 

The  case  where  ligand-receptor  associations  were  eliminated  shows  an  equilibrium 
value  of  975  complexes  reached  at  0.5  seconds.  This  equilibrium  is  much  larger  than 
for  the  default  case  because  there  are  no  i-1  or  j-l  complexes  present  in  the  system. 
All  of  the  receptors  first  form  the  i-j  heterodimer  and  then  bind  bulk  free  ligand  to 
produce  the  signaling  complex. 

The  case  where  preformed  heterodimers  are  eliminated  behaved  in  a  similar  fashion 
to  the  default  case.  The  complex  (i-j)-l  reached  an  equilibrium  value  of  320  complexes 
at  10  milliseconds  and  the  i-1  and  j-l  complexes  reached  equilibrium  values  of  680 


complexes  in  each  case  a‘  a  time  of  9  milliseconds.  The  number  of  i-1  and  j-1  formed 
is  slightly  higher  than  in  the  default  case  because  there  are  no  i  and  j  receptors  in  the 
form  of  the  preformed  heterodimer  i-j. 
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7.2.4  Paracrine  Growth  Case 

Figure  11  shows  the  case  for  a  system  that  allows  formation  of  ligand-receptor  com¬ 
plexes  through  the  paracrine  growth  mechanism  only.  The  paracrine  growth  only  case 
shows  behavior  that  differed  greatly  from  the  previously  described  cases.  The  default 
paracrine  growth  case  shows  a  large  number  of  preformed  heterodimers  present  before 
the  complex  i-j-1  was  produced  to  any  appreciable  amounts.  The  complexes  i-1  and 
j-1  are  not  present  to  any  great  extent  in  this  case. 

i 

The  case  where  ligand-receptor  associations  were  eliminated  shows  behavior  that 
more  closely  matches  the  two  previously  mentioned  cases.  The  complex  i-j-1  reached 
an  equilibrium  value  of  975  complexes  at  a  time  of  1  seconds  closely  matching  the 
equilibrium  number  and  time  of  the  autocrine  only  and  autocrine/ paracrine  cases. 
The  main  difference,  however,  is  the  presence  of  a  great  number  of  the  preformed  het¬ 
erodimer  at  60  milliseconds.  The  maximum  value  reached  is  just  over  900  complexes 
which  quickly  recedes  as  the  complex  i-j-1  is  produced. 

The  case  where  the  preformed  heterodimer  formation  is  eliminated  shows  un¬ 


usual  results  in  that  the  complex  i-j-1  is  produced  without  any  appreciable  numbers 
of  i-1  or  j-1  complexes  being'  present,  which  were  present  in  the  autocrine  and  au¬ 
tocrine/paracrine  systems. 
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7.2.5  Discussion  of  Results  for  the  Complex  Model  with  Identical  Recep¬ 
tors 

The  autocrine  and  autocrine/paracrine  systems  showed  almost  identical  behavior  and 
will  be  treated  together.  It  is  of  interest  to  note  that  in  both  of  these  systems,  the  case 
where  preformed  heterodimers  were  removed  from  the  model  most  closely  resembled 
the  default  data  case  where  all  models  were  present,  both  equilibrium  numbers  of 
complexes  formed  and  the  time  at  which  equilibrium  was  achieved.  For  these  systems 
then,  the  formation  of  the  signaling  complex  i-j-1  must  occur  via  the  affinity  conversion 
model.  The  case  where  ligand-surface  interactions  have  been  eliminated  shows  that 
indeed  the  signaling  complex  (i-j)-l  may  be  formed  by  the  preformed  heterodimer 
model.  However,  in  this  case  the  equilibrium  number  of  i-j-1  is  reached  at  a  later  time 
than  in  the  autocrine  and  autocrine/ paracrine  cases. 

This  allows  the  following  generalizations  to  be  made  about  the  nature  of  the  system 
when  applied  to  the  interleukin-2  problem.  The  complexes  i-1  and  j-1  must  be  formed 
faster  than  the  preformed  heterodimer.  This  is  shown  by  the  fact  that  the  equilibrium 
number  of  complexes  is  obtained  at  10  milliseconds  for  the  default  and  no  preformed 
heterodimer  cases  and  at  0.5  seconds  for  the  case  where  ligand-receptor  interactions 
have  been  eliminated.  The  preformed  heterodimer  model  is  suppressed  in  the  default 
case  because  the  i  and  j  receptors  bind  ligand  to  become  the  complexes  i-1  and  j-1 
much  faster  than  they  Hr  d  *'ach  other  to  become  the  preformed  heterodimer.  These 
complexes  then  diffurr .  >  .it  randomly  until  they  find  the  other  type  of  receptor  and 


form  the  signaling  complex  i-j-1. 

The  paracrine  growth  model  is  unique  in  that  no  non-specifically  bound  ligand  is 
initially  present  on  the  surface  of  the  cells.  The  default  case  shows  many  preformed 
heterodimers  before  the  production  of  the  signaling  complex  because  the  diffusion 
of  the  i  and  j  receptors  leads  only  to  the  formation  of  preformed  heterodimer.  The 
complexes  i-1  and  j-1  are  produced  only  by  the  ligand-surface  mechanism,  which  was 
shown  in  the  simple  model  to  lead  to  a  small  equilibrium  number  of  these  complexes 
when  compared  to  the  membrane- bound  mechanism  (Figure  3a).  The  results  of  the 
paracrine  only  case  reinforce  the  concept  that  the  membrane-bound  interactions  are 
the  predominant  interactions,  even  though  they  take  place  over  longer  time  periods 
than  ligand-surface  interactions.  Thus  the  autocrine/paracrine  model  shows  the  re¬ 
sults  obtained  for  the  autocrine  only  model,  where  ligand  surface  interactions  are  not 
present,  and  not  the  paracrine  only  model,  where  these  ligand  surface  interactions 
predominate. 

7.2.6  Comparison  of  Kinetic  Results  to  Actual  Experimental  Values 

Dissociation  rate  constants  were  measured  for  the  species  i-1,  g-1,  i-j,  and  (i-j)-l  and 
are  compared  to  experimentally  obtained  values  [14]  in  Table  1.  The  first-order  disso¬ 
ciation  plots  for  these  species  are  given  in  Figure  12.  The  association  rate  constant  for 
the  species  (i-j)-l  obtained  from  the  default  data  run  is  compared  to  an  experimentally 
obtained  value  [14]  in  Table  1. 


7.3  The  Complex  Model  with  Non-Identical  Receptors 


7.3.1  Experimental  Conditions 

All  runs  were  conducted  in  a  manner  similar  to  that  for  the  case  of  identical  receptors 
in  that  autocrine  only,  paracrine  only,  and  autocrine/paracrine  systems  were  exam¬ 
ined.  The  difference  was  that  the  probability  parameters  for  the  two  receptors  were 
changed  so  that  the  interleukin-2  ligand  bound  to  the  IL-2a  receptor  with  a  higher 
probability  than  it  bound  to  the  IL-2/?  receptor. 

7.3.2  Autocrine  Growth  Case 

Figure  13  shows  the  ligand-receptor  association  curves  for  the  autocrine  growth  only 

i 

case.  For  the  default  case,  it  is  of  interest  to  note  that  the  complex  i-1  is  formed  faster 
than  the  complex  j-1  and  that  both  of  these  species  reach  an  equilibrium  value  of  560 
complexes  at  times  of  1  and  5  milliseconds  respectively.  The  complex  (i-j)-l  reaches 
an  equilibrium  valu»,  of  380  complexes  at  a  time  of  10  milliseconds. 

The  case  where  ligand-receptor  associations  have  been  eliminated  shows  an  equi¬ 
librium  value  of  the  number  of  (i-j)-l  complexes  of  980  reached  at  a  time  of  0.5  seconds. 
This  is  because  no  i  and  j  receptors  are  used  up  in  the  form  of  the  complexes  i-1  or 
j-1,  thus  all  available  i  and  j  receptors  can  be  used  to  form  the  preformed  heterodimer 

>*j- 

The  case  where  preformed  heterodimers  have  been  eliminated  shows  behavior  that 
resembles  that  of  the  default  case.  The  complex  i-1  is  formed  faster  than  the  complex 


j-1  and  both  of  these  complexes  reach  equilibrium  at  1  milliseconds  and  5  milliseconds 
respectively.  The  numbers  reached  (625)  in  both  cases  is  greater  than  the  default 
case  because  no  i  and  j  receptors  are  in  the  form  of  the  preformed  heterodimer. 
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7.3.3  Autocrine  and  Paracrine  Growth  Case 

The  behavior  of  the  autocrine/paracrine  system  is  shown  in  the  ligand-receptor  asso¬ 
ciation  curves  given  in  Figure  14.  These  ligand-receptor  curves  closely  resemble  those 
obtained  in  the  autocrine  only  case. 

The  default  data  case  shows  that  the  i-1  complex  is  produced  faster  than  the  j-1 
complex.  Both  complexes  reach  an  equilibrium  value  of  560  complexes  at  times  of  1 
and  5  milliseconds,  respectively.  The  complex  (i-j)-l  reaches  an  equilibrium  value  of 
380  complexes  at  10  milliseconds. 

The  case  were  ligand-receptor  associations  have  been  eliminated  show  the  large 
number  of  (i-j)-l  complexes  characteristic  of  this  case,  with  an  equilibrium  number 
of  980  complexes  being  reached  at  0.5  seconds.  The  increased  number  of  complexes 
can  again  be  attributed  to  the  fact  that  no  i  or  j  receptors  are  used  up  as  either  the 
i-1  complex  or  the  j-1  complex.  Elimination  of  the  preformed  heterodimer  produces 
behavior  that  resembles  that  of  the  default  case.  The  production  of  the  complex  i-1 
occurs  faster  than  the  production  of  j-1.  Both  complexes  reach  an  equilibrium  value  of 
625  complexes  at  times  of  1  and  5  milliseconds,  respectively.  The  increased  number 
of  i-1  and  j-1  complexes  can  again  be  attributed  to  the  absence  of  the  preformed 
heterodimer,  providing  more  free  i  and  j  receptors. 
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7.3.4  Paracrine  Growth  Case 

The  ligand-receptor  association  curves  for  these  trials  are  shown  in  Figure  15.  The 
paracrine  growth  only  case  showed  behavior  that  differed  greatly  from  the  previously 
described  cases.  The  default  paracrine  growth  case  showed  a  large  number  of  pre¬ 
formed  heterodimers  present  before  the  complex  i-j-1  was  produced  to  any  appreciable 
amounts.  The  complexes  i-1  and  j-1  were  not  present  to  any  great  extent  in  this  case, 
except  at  long  times  where  dissociation  of  the  complex  (i-j)-l  gave  rise  to  these  species. 

The  case  where  ligand-receptor  associations  were  eliminated  shows  behavior  that 
more  closely  matches  the  two  previously  mentioned  cases.  The  complex  i-j-1  reached 
an  equilibrium  value  of  975  complexes  at  a  time  of  1  seconds  closely  matching  the  equi¬ 
librium  number  and  time  of  the  autocrine  only  and  autocrine/paracrine  cases.  The 
main  difference,  however,  the  presence  of  the  maximum  number  of  the  preformed 
heterodimer  at  60  milliseconds.  The  maximum  value  reached  was  just  over  900  com¬ 
plexes.  This  number  of  preformed  heterodimer  was  quickly  reduced  as  the  complex 
i-j-1  was  produced. 

The  case  where  the  preformed  heterodimer  formation  is  eliminated  shows  unusual 
results  in  that  the  complex  i-j-1  is  produced  without  any  appreciable  numbers  of  i-1 
or  j-1  complexes  being  present,  which  were  present  in  large  numbers  in  the  autocrine 
only  and  autocrine/paracrine  systems.  Here  the  complexes  i-1  and  j-1  are  only  present 
at  long  times  due  to  the  dissociation  of  (i-j)-l. 


7.3.5  Discussion  of  Results  for  the  Complex  Model  with  Non-  Identical 
Receptors 

The  results  for  the  complex  model  with  non-identical  receptors  closely  resemble  the 
results  obtained  for  the  complex  model  with  identical  receptors.  The  autocrine  and 
autocrine/paracrine  systems  both  behaved  very  similarly  and  will  be  examined  to¬ 
gether.  In  these  cases  the  affinity  conversion  model  is  again  preferred  over  the  pre¬ 
formed  heterodimer  models.  This  is  shown  by  the  fact  that  the  system  behaved  the 
same  in  both  the  default  case  and  the  case  where  preformed  heterodimers  had  been 
eliminated.  In  these  cases  the  complex  (i-j)-l  was  formed  by  10  milliseconds  compared 
to  0.5  seconds  for  the  case  where  ligand-receptor  associations  had  been  eliminated. 
The  formation  of  the  complexes  i-1  and  j-1  is  more  rapid  than  the  formation  of  the  pre¬ 
formed  heterodimer  and  thus  the  complex  i-j-1  is  produced  by  the  affinity  conversion 
mechanism  before  it  can  be  produced  by  the  preformed  heterodimer  mechanism. 

As  was  the  case  for  the  complex  model  with  non-identical  receptors,  the  paracrine 
growth  model  is  unique  in  that  no  non-specifically  bound  ligand  is  initially  present  on 
the  surface  of  the  cells.  The  default  case  shows  many  preformed  heterodimers  before 
the  production  of  the  signaling  complex  because  the  diffusion  of  the  i  and  j  recep¬ 
tors  leads  only  to  the  formation  of  preformed  heterodimer.  A  small  number  of  the 
complexes  i-1  and  j-1  are  present  because  these  complexes  are  formed  ligand-surface 
associations  only.  As  shown  in  the  simple  model,  this  leads  to  a  small  equilibrium 
number  of  these  complexes  when  compared  to  the  membrane-bound  mechanism  (Fig- 
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ure  3a). 

The  autocrine/paracrine  model  shows  the  results  obtained  for  the  autocrine  only 
model,  where  ligand  surface  interactions  are  not  present,  and  not  the  paracrine  only 
model,  where  these  ligand  surface  interactions  predominate.  This  again  shows  that 
the  membrane-bound  diffusion  interactions  are  the  most  significant  interactions  at 
long  times. 

8  Conclusion 

The  effectiveness  of  the  application  of  a  diffusion-kinetics  system  to  ligand-receptor 
systems  according  to  the  method  of  Waite  [1]  has  been  demonstrated.  The  simple 

i 

model  has  been  used  to  describe  the  nature  of  these  ligand-receptor  interactions  and 
prove  the  validity  of  theoretical  results  through  comparison  with  actual  experimen¬ 
tally  obtained  values.  The  complex  model  has  been  used  to  show  that  the  affinity 
conversion  model  is  the  preferred  mechanism  for  the  formation  of  the  interleukin-2  sig¬ 
naling  complex  according  to  theoretical  diffusion  kinetics.  The  theoretically  obtained 
rate  constants  and  half-lives  given  in  Table  3  differ  somewhat  from  experimentally 
obtained  values  [14]  indicating  that  better  estimates  of  parameters  such  as  diffu¬ 
sion  coefficients  and  radii  are  needed  before  the  model  simulates  the  behavior  of  the 
interleukin-2  system  with  complete  accuracy. 


FIGURE  LEGENDS 


FIGURE  1.  Ligand  -  receptor  association  curves.  The  association  curves  for  the 
IL-2-a  receptor  complex  (shaded  box  for  r  =  lxlO“6m  and  +  for  r  =  lxlO“5m)  and 
non-specifically  bound  ligand  (*  for  r  =  lxlO“6m  and  empty  box  for  r  =  lxlO“5m): 
A.  The  association  curves  used  to  determine  the  association  rate  constants  (listed  in 
Table  3)  for  the  IL-2-a  receptor  complex  (shaded  box  for  r  =  lxlO~6m  and  +  for  r  = 
lxlO-5m):  B.  The  association  curves  for  the  binding  of  IL-2  by  the  a  receptor  showing 
switch  from  the  predominance  of  bulk  to  surface  associations  to  the  predominance 
of  membrane-bound  diffusional  associations  at  3.0  msec:  C.  The  association  rate 
constants  for  the  binding  of  IL-2  by  the  a  receptor  are  determined  from  the  initial 
slope  as  shown  by  Landgraf,  et  al.  [14]. 

FIGURE  2.  Ligand  -  receptor  association  curves.  The  association  of  IL-2  and  a 
receptors  for  default  data  (shaded  box),  Pi-i  =  0  (+),  and  P2,_j  =  0  (*):  A.  The 
association  curves  used  to  determine  the  association  rate  constants  (listed  in  Table 
3)  for  the  IL-2-a  receptor  complex  for  default  data  (shaded  box),  P,_j  =  0  (+), 
and  P2l_/  =  0  (*):  B.  The  association  curves  used  to  determine  the  association  rate 
constants  (listed  in  Table  )  for  the  IL-2-a  receptor  complex  for  default  data  (shaded 
box),  P,_;  =  0  (+),  and  P2,_/  =  0  (*):  C.  The  association  rate  constants  for  the 


binding  of  IL-2  by  the  a  receptor  are  determined  from  the  initial  slope  as  shown  by 
Landgraf,  et  al.[14]. 

FIGURE  3.  Ligand  -  receptor  association  curves.  The  association  of  IL-2  and  a 
receptors  for  Pi-i  =  0  and  r  =  lxlO“6m  (shaded  box),  Pn-i  =  0  and  r  =  lxlO~6m 
(+),  Pi-t  —  0  and  r  =  lxl0“8m  (*),  and  PM  =  0  and  r  =  lxlO“8m  (empty  box): 
A.  The  association  curves  used  to  determine  the  association  rate  constants  (listed  in 
Table  )  for  Pi-i  =  0  and  r  =  lxlO“6m  (shaded  box)  and  P,_j  =  0  and  r  =  lxlO“8m 
(+):  B.  The  association  curves  used  to  determine  the  association  rate  constants 
(listed  in  Table  3)  for  P2,_/  =  0  and  r  =  lxl O'® m  (shaded  box)  and  Pj,_/  =  0  and 
r  =  lxlO~8m  (+):  C.  The  association  rate  constants  for  the  binding  of  IL-2  by  the  a 
receptor  are  determined  from  the  initial  slope  as  shown  by  Landgraf,  et  al.  [14]. 

FIGURE  4.  Ligand  •  receptor  association  curves.  The  relationship  between  the 
number  of  IL-2-a  receptor  complexes  and  the  number  of  non-specifically  bound  ligand 
is  shown  for  a  cell  radius  of  lxlO”6m  (shaded  box  for  IL-2-a  and  +  for  g-1):  A.  The 
relationship  between  the  number  of  IL-2-o  receptor  complexes  and  the  number  of 
non-specifically  bound  ligand  is  shown  for  a  cell  radius  of  lxlO~8m  (shaded  box  for 
IL-2-or  and  +  for  g-1):  B. 


FIGURE  5.  Ligand  •  receptor  association  and  dissociation  curves.  The  association 
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curves  for  the  binding  of  IL-2  by  the  a  receptor  for  Pg-i  =  0,  P,_/  =  0,  Pj-i  =  0,  and 
r  =  lxl0“6m  (shaded  box)  and  P,_j  =  0,  Pa,_/  =  0,  P^j-i  —  0,  and  r  =  lxl0-6m 
(+):  A.  The  association  curves  for  the  binding  of  IL-2  by  the  a  receptor  for  P#_/  =  0, 
Pi-i  =  0,  Pj-i  =  0,  and  r  =  lxl0_5m  (shaded  box)  and  Ps-i  =  0,P2,_/  =  0,  Pij-i  =  0, 
and  r  =  lxl0“5m  (+):  B. 

FIGURE  6.  Kinetics  of  ligand  dissociation.  The  dissociation  of  IL-2  at  non-specific 
sites:  shaded  box,  at  the  a  receptor  on  cells  of  radius  lxl0-6m:  empty  box,  at  the 
a  receptor  on  cells  of  radius  lxl0"sm:  vertical  line.  The  dissociation  rate  constants 
determined  from  the  slopes  of  the  curves  are  listed  in  Table  3  and  the  half-lives 
determined  from  the  dissociation  rate  constants  are  listed  in  Table  3. 

FIGURE  7.  Ligand  -  receptor  association  curves.  The  association  curves  for  the 
binding  of  IL-2  by  the  a  receptor  for  r  =  1  x  1 0_6mandZ?j^a„d  =  lxlO~10m2s_1  (shaded 
box),  r  =  lxl0~6m  and  Dusan<i  —  lxl0-15mJs-1(+),  and  r  =  lxl0~6m  and  Dugand  — 
lxl0_5m2s~1(*):  A.  The  association  curves  for  the  binding  of  IL-2  by  the  a  recep¬ 
tor  for  r  =  1x10 ~sm  and  Dugand  =  lxl0_10m2s”1(shaded  box),  r  =  lxl0~5m  and 
Dugand  —  lxl0“15m2s“1(+),  and  r  =  lxl0"5m  and  Dngand  —  lxl0-5m2s_I(*):  B. 

FIGURE  8.  Ligand  -  receptor  association  curves.  The  association  curves  for  the 
binding  of  IL-2  by  the  a  receptor  for  [ligand\0  —  0  and  number  g-l0  =  100000(shaded 


box)  and  for  [ligand\0  =  0,  number  g  -  l0—  100000,  and  P,_j  =  0  (+). 

FIGURE  9.  Ligand  •  receptor  association  curves  for  the  autocrine  growth  system 
and  identical  receptors.  The  legend  for  all  curves  is  shaded  box:  i-1,  +:  i-j,  *:(i-j)-l, 
empty  box:  (i-j)...l.  Default  Data:  A.  Elimination  of  ligand-receptor  associations:  B. 
Elimination  of  preformed  heterodimer:  C. 

FIGURE  10.  Ligand  -  receptor  association  curves  for  the  autocrine/paracrine  growth 
system  and  identical  receptors.  The  legend  for  all  curves  is  shaded  box:  i-1,  -f :  i-j,  *:(i- 
j)-l,  empty  box:  (i-j)--.I.  Default  Data:  A.  Elimination  of  ligand-receptor  associations: 
B.  Elimination  of  preformed  heterodimer:  C. 

FIGURE  11.  Ligand  -  receptor  association  curves  for  the  paracrine  growth  system 
and  identical  receptors.  The  legend  for  all  curves  is  shaded  box:  i-1,  +:  i-j,  *:(i-j)-l, 
empty  box:  (i-j)...l.  Default  Data:  A.  Elimination  of  ligand-receptor  associations:  B. 
Elimination  of  preformed  heterodimer:  C. 

FIGURE  12.  First  order  kinetic  plots  for  the  dissociation  of  g-1  (shaded  box),  i-j  (+), 
i-1  (empty  box),  and  (i-j)-l  (*).  The  dissociation  constants  are  listed  in  Table  3. 

FIGURE  13.  Ligand  -  receptor  association  curves  for  the  autocrine  growth  system 
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and  non-identical  receptors.  The  legend  for  all  curves  is  shaded  box:  i-1,  +:  j-1, 
*:i-j,  empty  box:  (i-j)-l,  x:(i-  j)...l.  Default  Data:  A.  Elimination  of  ligand-receptor 
associations:  B.  Elimination  of  preformed  heterodimer:  C. 

FIGURE  14.  Ligand  -  receptor  association  curves  for  the  autocrine/paracrine  growth 
system  and  non-identical  receptors.  The  legend  for  all  curves  is  shaded  box:  i-1,  +:  j-1, 
*:i-j,  empty  box:  (i-j)-l,  x:(i-j)...l.  Default  Data:  A.  Elimination  of  ligand-  receptor 
associations:  B.  Elimination  of  preformed  heterodimer:  C. 

FIGURE  15.  Ligand  -  receptor  association  curves  for  the  paracrine  growth  system 
and  non-identical  receptors.  The  legend  for  all  curves  is  shaded  box:  i-1,  +:  j-1, 
*:i-j,  empty  box:  (i-j)-l,  x:(i-  j)...l.  Default  Data:  A.  Elimination  of  ligand-receptor 
associations:  B.  Elimination  of  preformed  heterodimer:  C. 


References 


[1]  Waite,  B.A.  1993.  A  Diffusion-Kinetics  Model  for  Receptor/Ligand  Interactions 
at  Cell  Surfaces.  Manuscript. 

[2]  Waite,  B.A.,  and  J  D.  Stewart.  1992.  An  idealized  dynamical  model  of  simple 
diffusional  interactions  between  macromolecules  and  between  macromolecules 
and  surfaces.  Mathematical  Biosciences,  in  press. 

[3]  Darnell,  James;  Lodish,  Harv  ;;  Baltimore,  David,  Molecular  Cell  Biology ;  Sci¬ 
entific  American  Books:  New  York,  1986;  p.  667-  669. 

[4]  Singer,  S.J.  and  Micolson,  G.L.,  "The  Fluid  Mosaic  Model  of  the  Structure  of 
Cell  Membranes,”  Science  173,  p.  720-731  (1972). 

[5]  Golub,  Edward  S.  and  Green,  Douglas  R,  Immunology:  A  Synthesis,  2nd  Edition; 
Sinauer  Associates,  Inc:  Sunderland,  Massachusetts,  1991,  p.  456. 

[6]  Smith,  K.A.,  Science  240,  p.1169-1176  (1988). 

[7]  Caligiuri,  M.,  Zmuidzinas,  A.,  Manley,  T.J.,  Levine,  H.,  Smith,  K.A.,  and  Ritz, 
J.,  J.  Exp.  Med.  171,  p.  1509-1526  (1990). 

[8]  Golub,  Edward  S.  and  Green,  Douglas  R,  Immunology:  A  Synthesis,  2nd  Edition; 
Sinauer  Associates,  Inc:  Sunderland,  Massachusetts,  1991,  p.  453. 


52 


[9]  Smith,  K.A.,  Ann.  Rev.  Cell  Biol.  5,  p.  397-425  (1989). 

[10]  Golub,  Edward  S.  and  Green,  Douglas  R,  Immunology:  A  Synthesis,  2nd  Edition; 
Sinauer  Associates,  Inc:  Sunderland,  Massachusetts,  1991,  p.  455. 

[11]  Golub,  Edward  S.  and  Green,  Douglas  R,  Immunology:  A  Synthesis,  2nd  Edition; 
Sinauer  Associates,  Inc:  Sunderland,  Massachusetts,  1991,  p.  454. 

[12]  Wang,  H.M.,  and  Smith,  K.A.,  J.  Exp.  Med.  166,  p.  239-243  (1989). 

[13]  Teshigawara,  K.,  Wang,  H.M.,  Kato,  K.,  and  Smith,  K.A.,  J.  Exp.  Med.  165  p. 
223-238  (1987). 

[14]  Landgraf,  B.E.,  D.P.  Williams,  J.R.  Murphy,  T.R.  Sana,  K.A.  Smith,  and  T.L. 
Ciardelli.  1992.  Working  Manuscript. 

[15]  Goldstein,  B.,  D.  Jones,  I.G.  Kevrekidis,  and  A.S.  Perelson.1992.  Evidence  for 
p55-p75  heterodimers  in  the  absence  of  interleukin-2  from  Scatchard  Plot  Anal¬ 
ysis.  International  Immunology,  4:23-32. 

[16]  Kamio,  M.,  T.  Uchiyama,  N.  Arima,  K.  Itoh,  T.  Ishikawa,  T.  Hori,  and  H. 
Uchino.  1990.  Role  of  chain-IL-2  complex  in  the  formation  of  the  ternary  complex 
of  IL-2  and  high-affinity  IL-2  receptor.  International  Immunology,  2:521-530. 

[17]  Saragovi,  H.  and  T.R.  Malek.  1988.  Direct  identification  of  the  murine  IL-2 
receptor  p55-p75  heterodimer  in  the  absence  of  IL-2.  Journal  of  Immunology, 


141:476-482 


52 


a 


[18]  Grant,  A.J.,  E.  Roessler,  G.  Ju,  M.  Tsudo,  K.  Sugamura,  T.A.  Waldmann.  1992. 
The  interleukin-2  receptor  (IL-2R):  The  IL-2R  subunit  alters  the  function  of 
the  IL-2R  subunit  to  enhance  IL-2  binding  and  signaling  by  mechanisms  that 
do  not  require  binding  of  IL-2  to  IL-2R  subunit.  Proc.  Natl.  Acad.  Sci.,  USA, 


89:2165-2169 


53 


Table  1 


species 

number 

1 

2 

3 

4 

5 

6 
7 
i 

9 

10 
11 


Table  1A 

Species  of  the  Simple  Model 

species  description 

symbol 

I  bulk  free  intcrlcukJn-2  ligand 

i  11-2  alpha  receptor 

j  11-2  beta  receptor 

i-1  ligand  bound  to  ]|-2  alpha  receptor 

j-1  ligand  bound  to  11-2  beta  receptor 

g-1  ligand  non-specifically  bound  to  cell  surface 

i~l  ligand  rebounding  from  the  11-2  alpha  receptor 

j— 1  ligand  rebounding  from  the  11-2  beta  receptor 

J— 1  ligand  rebounding  from  a  non-specific  binding  site 

»— g-1  non-spccifcally  bound  ligand  rebounding  from  11-2  alpha  receptor  on  the  cell  membrane 
j— g-1  non-spccifcally  bound  ligand  rebounding  from  11-2  beta  receptor  on  the  cell  membrane 


Table  IB 

Species  of  the  Complex  Model 

The  complex  model  includes  all  of  the  species  of  the  simple  model  and  the  following  additional  species 


species 

species 

number 

symbol 

12 

‘-j 

13 

i—j 

14 

Uj-I 

15 

j — i-l 

16 

(i-j)-l 

17 

(i-j)-l 

description 

heterodimer  of  11-2  alpha  and  11-2  beta  receptors 

rebounding  encounter  between  11-2  alpha  and  11-2  beta  receptors 

11-2  alpha  receptor  rebounding  from  the  11-2  beta  receptor-ligand  complex 

11-2  beta  receptor  rebounding  from  the  11-2  alpha  receptor-ligand  complex 

intcrieukin-2/11-2  alpha  receptor/11-2  beta  receptor  signalling  complex 

ligand  rebounding  from  the  11-2  alpha  receptor/11-2  beta  receptor  betcrodimcr 


T 


Table  2  Experimental  Parameter* 

Valne 

Complex  Model 

Simple 

Idcitieil 

Noa-ldeat 

Parameter 

Model 

Recepton 

Recepton  Dcccriptioa 

pa  . 

a* 

03 

03  Probability  o£  binding  tor  i+l  via  ligaad-aaifacc  aeckaaiam 

Pjl 

a» 

03 

01  Probability  of  binding  for  j+i  via  ligaod-aaifaeeaccbaaiaa 

P*l 

as 

as 

OS  Probability  cl  btadiag  for  g+l  via  ligaed-anrfaec  mechaaiaa 

P2il 

as 

as 

OS  Probe bflity  of  binding  for  i + 1  via  ncmbnac-booad  dif tacioa  aeckaaiam 

n? 

as 

as 

02  Probability  of  binding  forj+l  via  mcnbnnc-boaad  diffuion  meckaaiam 

Pij 

a/a 

as 

OS  Probability  of  binding  for  i+ j  via  membnae-bonad  diffntioa  methaaiaa 

Pijl 

n/» 

as 

OS  Probability  of  binding  to  fora  (i-j)-l  via  a  ligand  larfacc  or  acabraae-booad 
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Title  3:  Summary  of  Kinetic  Dito 


Association  Rite  Constants 
l/(sccond*M) 


radius 

Simple 

Model 

1  micron 

10  mieron 

1  mieron 

i-l 

i-l 

i-i 

default 

4.52E+06 

4J2E+07 

4J2E+06 

pil=0 

4.42E+06 

^lE+07 

4.54E+06 

p2il-0 

2.47E+05 

L87E+04 

&56E+06 

Dissociation  Rate  Constants 

l/seeonds 

Simple 

i-l 

2.02E-03 

8-1 

1.00E-03 

Model 

Complex 

2.02E-03 

9J9 OE-04 

Model 

Landgraf 

430E-02 

n/a 

Previously  • 

Complex 

Prev. 

Reported[14] 

Model 

Reported  [14J 

• 

1  micron 

i-l 

(i-J>l 

1.20B+07 

•  &79E+10 

1.2QE+08 

n/a 

n/a 

n/a 

n/a 

n/a 

n/a 

i-j  (i-j>l 

n/a  n/a 


&10E-04  4.04E-03 


n/a  1.60E-04 


Half-Lives 

seconds 

i-1  |-1  i-j  (i-j)-l 

Simple  3.43E+02  6JT2E+02  n/a  n/a 

Model 

Complex  3.43E+02  &98E+02  8^6E+02  1.71E+02 

Model 

Landjraf  1.54E+01  n/a  n/a  4J3E+Q3 


10  IE-08  IE-06  IE-04  IE-02  1E+00  1E+02  1E+04 

IE-09  IE-07  IE-05  IE-03  IE-01  1E+01  1E+03 

Seconds 


Number 


Number  of  Complexes 


72 


Number  of  Complexes  Present 


Number  of  Complexes  Present 


75 
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Number  of  Complexes  Present 


Figure  IOC 
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I 


Number  of  Complexes  Present 


Figure  13A 
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E-06  IE-04  IE-02  1E+00 

IE-05  IE-03  IE-01  1E-» 

Seconds 
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E-09  IE-07  IE-05  IE-03  IE-01  1E+01  1E+03 

Seconds 


I 


IE-08  IE-06  IE-04  IE-02  1E+00  1E+G2  1E+04 

E-09  IE-07  IE-05  IE-03  IE-01  1E+01  1E+03 

Seconds 


